# 此程序可将CGCS2000坐标系下的平面坐标值转换为此坐标系下的经纬度
from pyproj import Proj, transform

# CGCS2000 坐标数据示例
cgcs2000_coordinates = [
    (773118.408,3998395.095) # x, y 坐标 示例数据为安阳地图录制文件第一行
    # 你可以继续添加更多的坐标
]
 
# 定义 CGCS2000 坐标系
cgcs2000 = Proj(proj='lcc', lat_1=30, lat_2=40, lat_0=35, lon_0=105, x_0=0, y_0=0, 
                 ellps='GRS80', towgs84='0,0,0,0,0,0,0', no_defs=True)
 
# 定义 WGS 84 经纬度坐标系
wgs84 = Proj(proj='longlat', datum='WGS84')
 
def convert_coordinates(cgcs2000_x, cgcs2000_y):
    """
    将 CGCS2000 坐标转换为经纬度
    :param cgcs2000_x: CGCS2000 坐标 x 分量
    :param cgcs2000_y: CGCS2000 坐标 y 分量
    :return: 经纬度元组 (经度, 纬度)
    """
    longitude, latitude = transform(cgcs2000, wgs84, cgcs2000_x, cgcs2000_y)
    return longitude, latitude

# 测试转换功能
for x, y in cgcs2000_coordinates:
    longitude, latitude = convert_coordinates(x, y)
    print(f"CGCS2000 坐标: (X: {x}, Y: {y}) -> 经度: {longitude}, 纬度: {latitude}")